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Monolayer Topological Insulators: Silicene, Germanene and Stanene 
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We report the recent progress on the theoretical aspects of monolayer topological insulators including silicene, 
germanene and stanene, which are monolayer honeycomb structures of silicon, germanium and tin, respectively. 
They show quantum spin Hall effects in nature due to the spin orbit interaction. The band gap can be tuned by 
applying perpendicular electric field, which induces a topological phase transition. We also analyze the topolog¬ 
ical properties of generic honeycomb systems together with the classification of topological insulators. Phase 
diagram of topological insulators and superconductors in honeycomb systems are explicitly determined. We also 
investigate topological electronics including a topological field-effect transistor, the topological Kirchhoff’s law 
and the topological spin-valleytronics. 


I. INTRODUCTION 


Monolayer materials are one of the most active fields of 
condensed matter physics. Graphene, monolayer honeycomb 
structure of carbon atoms, is the first experimentally real¬ 
ized monolayer material. Its low-energy band structure is de¬ 
scribed by the Dirac theory, which results in various novel 
physical propertieJ®^. The success of graphene evokes an 
extensive search for other monolayer materials. In particu¬ 
lar, monolayer topological materials are fascinating, realizing 
topological insulators and topological superconductors. 

A natural question is whether other monolayer honeycomb 
systems purely made of one kind of atoms are possible. It is 
shown that monolayer honeycomb systems made of silicon, 
germanium and tin are possible, which are named silicene, 
germanene and stanene, respectively. Silicene is named af¬ 
ter the combination of silicon and suffix "ene", which means 
the sp 2 bonding structurJ^. Germanene is also named after 
germanium -b "ene". Stanene is named after the Latin word 
"stannum" for tirP. 

Silicene, Germanene and stanene are expected to be topo¬ 
logical insulators. Topological insulator (TI) is a distinctive 
state of matter indexed by topological numbers, and charac¬ 
terized by an insulating gap in the bulk accompanied by topo¬ 
logically protected gapless edgeJ^. Thus the physics of these 
materials is located at the confluence of graphene and topo¬ 
logical insulators, which results in very rich physics. 

Silicene and ge rmanene are proposed by first-principles 
calculationP^^^®21 where their stability and the emergence of 
the Dirac cone are predicted by first-principles calculations. 
It is shown that the band gap of silicene is electrically con¬ 
trollable by applying perpendicular el ectric field to silicene 
plane by the first-principle calculation^I^ES and based on the 
Dirac theor}P^. A first-principles calculation shows that the 
Dirac cones are hidden in silicene fabricated on the Ag sub¬ 
strate due to the strong hybridization between silicene and Ag 
substrat^^. 

First suggestive observationof silicene were reported 
in 2010. Silicene has been grown on various substrates 
such as the Ag substratJ^^M^ and the ZrB 2 substratJ^, the Ir 
substratJ^ and the M 0 S 2 substratJ^. Germanene is synthe¬ 
sized on the Au substratand Pt substratThere are sev¬ 
eral experiments on silicenJ^^®^. In particular, silicene was 
demonstrated in 2015 to act as a field-effect transistor at room 



FIG. 1: (a) The honeycomb structure, made of two fundamental vec¬ 
tors ai and a2, consists of two sublattices made of A and B sites. A 
dotted rectangular represents a unit cell, (b) The reciplocal lattice is 
also a honeycomb lattice. A dotted rectangular represents a unit cell, 
which contains two inequivalent points K and K'. 


temperaturJ^. There are reviews on experimental aspects of 
silicenJ^^Jf^. 

There are several proposal on realizing free-standing like 
silicene on substrates by first-principles calculations. Sil¬ 
icene will be grown on graphen^, hexagonal boron-nitridJ^, 
hydrogen-processed Si(lll) surfacJ^^, hydrogen-processed 
Ge(l 11) surfacJ^, Cl-passivated Si(l 11) and clean CaF 2 (l 11) 
surfacesolid argorF^I , between bilayer graphen^ and in- 
tercalating alkali metal atoms between silicene and the metal 
substrate^. 


II. GRAPHENE AND SILICENE 


The basic structure of graphene and silicene is a honeycomb 
lattice generated by the fundamental translational vectors ai 
and a 2 . It consists of two triangular sublattices made of in¬ 
equivalent lattice sites A and B [Figj^a)]. The reciprocal 
lattice is also a honeycomb lattice in the momentum space 
[Figj^b)], which constitutes the Brillouin zone. 


A. Graphene 

Graphene is described by the simplest tight-binding model 
on a honeycomb lattice [Fig|^, 

Hq = —t c\^Cjs^ (2.1) 

{hj)s 
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FIG. 2: Band structure of graphene, (a) Six valleys are seen in this 
figure, (b) The gap is closed at the K and K' points, where the 
band structure looks like a cone. It is called the Dirac cone because 
the dispersion is linear. Among six cones only two are physically 
inequivalent. 


where c\^ creates an electron with spin polarization s =ti at 
site i, {i,j) runs over all the nearest neighbor hopping sites, 
and t is the transfer energy. By diagonalizing the Hamiltonian 
we obtain the band structure, which we illustrate in Fig|^a). 
It consists of valleys or cones near the Fermi surface. The 
cones touch the Fermi surface at two inequivalent points, that 
is, the K and K' points in the Brillouin zone [Fig|^b)]. 

We are interested in physics near the Fermi energy. To de¬ 
rive the relevant Hamiltonian, we rewrite as 


in the momentum space, with 

/ (fe) = cos —. 


CBs ) 
( 2 . 2 ) 


(2.3) 


The energy spectrum is obtained as 


_ , / akr ySak.j ^ ak^ 

E (k) = ty 1 + 4 cos cos ——^ + 4 cos^ . (2.4) 

^ z z z 

The gap closes at the point defined by 


1 / Ttt 
^v = - ( 


with 


r] = ±. (2.5) 


The and K_ points are identical to the K and K' points, 
respectively. Because the dispersion relation is linear for ki 
0, they are also called the Dirac points: See Figsand [^a2). 

In the vicinity of the point, the Hamiltonian is approxi¬ 
mated by 

= E / ) , (2.6) 

with 

HJ = ivKr. + k,r,) = „ •*. ) 

(2.7) 

where r = {rx^Ty^Tz) is the Pauli matrix of the sublattice 

pseudospin for the A and B sites, and up = is the Fermi 
velocity with a being the lattice constant. The dispersion re¬ 
lation is linear for ki c:z0. We refer to as the Dirac Hamil¬ 
tonian at the Dirac point Ky. 



FIG. 3: (al) The lattice structure of graphene is planar, but (bl) that 
of silicene is buckled. Red and blue balls represent A and B sites. 
(a2) The band gap of graphene is closed, where the dispersion is lin¬ 
ear near the Fermi energy. (b2) The band gap of silicene is open. 
(a3) A fiat line connecting the K and K' points represents gapless 
flat edge modes in a nanoribbon. It contains 4-fold degenerate edge 
states for up/down-spin and left/right movers. (b3) Two lines con¬ 
necting the tips of the Dirac cones represents gapless edge modes of 
a nanoribbon. Each line contains 2-fold degenerate edge states. 


B. Silicene and tunable band gap 

The basic nature of silicene is described also by the tight- 
binding model There are two additional features mak¬ 

ing silicene essentially different from graphene. One is the 
presence of the spin-orbit interaction, which makes silicene a 
topological insulatoi®. The other is its buckled structure with 
a layer separation between the two sublattices [Fig|^bl)]. 
This freedom allows us to tune the gap by introducing a po¬ 
tential difference between the two sublattice^^^^. When we 
apply electric field perpendicular to silicene, the tight- 
binding Hamiltonian reads 

^ E + E 

( 2 . 8 ) 

is 

where ((i, j)) run over all the next-nearest neighbor hopping 
sites. The spin index stands for s for indices and for 
s = E within equations. It describes germanene and stanene 
as well. 

We explain each term, (i) The first term represents the usual 
nearest-neighbor hopping with the transfer energy t. (ii) The 
second term represents the effective SO coupling with Aso, 
where Uij = -fl if the next-nearest-neighboring hopping is 
anticlockwise and Uij = — 1 if it is clockwise with respect to 
the positive 2 ; axi^. (iii) The third term represents the stag¬ 
gered sublattice potential with jzi = -j-l (—1) for the A (B) 
sitJf^. Explicit values of these parameters are summarized 
in the Table [T| By diagonalizing the Hamiltonian by setting 
Ez = 0, we obtain the band structure illustrate as in Fig[^b2). 
The prominent feature is that the gap is open due to the SO 
interaction, and hence silicene is an insulator. A large SO in¬ 
teraction with Aso = 0.3eV is materialized in functionalized 
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i(eV) 

V 

a(A) 

Aso 

Ar 

£ 

0 

Graphene 

2.8 

9.8 

2.46 

10“^ 

0 

0 

90 

Silicene 

1.6 

5.5 

3.86 

3.9 

0.7 

0.23 

101.7 

Germanene 

1.3 

4.6 

4.02 

43 

10.7 

0.33 

106.5 

Stanene 

1.3 

4.9 

4.70 

43 

9.5 

0.33 

107.1 


TABLE I: The parameters charactering graphene, silicene and ger- 
manene . Here, vy is in the unit of 10^m/s, and Aso in the unit of 
meV. Ar is the Rashba SO interaction strength in the unit of meV: 
See LDf, i is the buckle height, while 0 is the bond angle. Taken 
from RefP. 


(a) trivial insulator 



(b) topological insulator 




FIG. 5: Berry curvature in (a) the trivial insulator and (b) the topo¬ 
logical insulator. It is strictly localized in the vicinity of the K and 
K' points. The Chern number, the integration of the Berry curva¬ 
ture over the Brilloin zone, is zero for a trivial insulator, while it is 
non-zero for a topological insulator. 


FIG. 4: Electrically tunable band gap and topological phase transi¬ 
tion of silicene. Silicene is a QSH insulator without electric field. By 
applying electric field, the band gap reduces and closes at the critical 
electric field ibF^cr. Above the critical electric field, silicene becomes 
a QVH insulator. 


cone is indexed by two parameters r] = ± and s = ±, the 
most general Dirac mass must have the following expression, 

= rjsXso — Ay + tjXh + sXsx, (2.13) 


stanen^, which will be a topological insulator at room tem¬ 
perature. 

The low-energy physic near the Fermi energy is described 
by the Dirac theory, which is constructed just as in the case 
of graphene. We rewrite the Hamiltonian ( |2.8| ) in the form of 
The Dirac Hamiltonian is explicitly given by 


H^ = 
where 

s 

with 


<1 

hupirjkx - iky) \ 

(2.9) 

hvF{rikx + iky) 

<1 

1 

7]sXso - ^Ez = 

-£{Ez - vsEa), 

(2.10) 

Ect = 

^so/£- 

(2.11) 


Note that A^ acts as the Dirac mass. The energy spectrum 
reads 

E (k) = ±yj {hv^kf + . (2.12) 

The gap is given by 2|A2| = 2i\Ez — ^sF^crl- 

It is important that the band gap is tunable by controlling 
external electric field Ez. The gap is open when Ez = ^. 
As \Ez\ increases, the gap become narrower [FigQ, and it 
closes at Ez = r]sEcv, where silicene is semimetallic just as 
in graphene. As \Ez \ increases further, the gap opens again. 


so that it has four independent parameters, Aso, Ay, Xh and 
Xsx^ We have already discussed the first two terms represent¬ 
ing the SO interaction and the sublattice staggered potential 
with Ay = ^Ez. We may write down the tight-binding terms 
that yield the fourth and fifth termJ^, 


. Xh 


{{id))s 


As'X ^ S^ic\^Cis. 
is 


(2.14) 


The fourth term describes the Haldane interaction induced by 
the photo-irradiation, where Xq = with Q the fre¬ 

quency and A the dimensionless intensit})^^®^. The fifth term 
describes the antiferromagnetic exchange magnetizatioiJ^. 

Here we note that there are a variety of 2D materials whose 
low-energy physics is described by the Dirac Hamiltonian 
( |2.9| ) with the Dirac mass ( 2.13| ). We call them general hon¬ 
eycomb systems. Examples are monolayer antiferromagnetic 
manganese chalcogenophosphates (MnPXs, X = S, Sej^Hand 
perovskite G-t^e antiferromagnetic insulators grown along 
[111] directiorPS. 

In what follows we analyze the Dirac Hamiltonian (2.9) 
with the Dirac mass ( |2.13| ). It can be positive, negative or 
zero. The band gap is given by 2| A^ |. 


III. TOPOLOGICAL PHASE TRANSITION 
A. Chern numbers 


C. Generalized Dirac mass terms 

There are actually other ways to control the band gap by 
introducing other interactions to silicene. Since each Dirac 


For any insulating state |t/; (k)) we may define a "gauge 
potential" in the momentum space by 

ak (fe) = -i (V’ (fc)l dk IV’ (fc)), 


(3.1) 
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which is properly called the Berry connection. Then we may 
define the "magnetic field" in the momentum spac, which is 
properly called the Berry curvature F (k), 

f) f) 

The Chern number is the integral of the Berry curvature F (k) 
over the first Brillouin zone, which is the total "magnetic fiux", 



(3.3) 


We have calculated the Berry curvature with the use of the 
tight-binding Hamiltonian ( |2.8| ), which we illustrate in Fig|^ 
The Berry curvature is strictly localized at the K and K' 
points. This feature remains unchanged even if we include 
the extra terms ( [2.14| ). Consequently, the Dirac Hamiltonian 
is valid to make a topological analysis to each valley, which 
is indexed by the spin s and the valley index r] = ±. 
Namely, it is possible to assign the Chern number to each 
valley. 

When the Hamiltonian is given by ( |2.9| ), the Berry curvature 
is explicitly calculated for each valley as 


F; (fe) = -r,- 




({hv^kf + (A^f^ 


3/2 ■ 


The Chern number is obtained as 


Cl = -/sgn(A^), 


(3.4) 


(3.5) 


where the Dirac mass is given by ( |2.13[ ). 

The Chern number is quantized as CJ ± |. It is insensi¬ 

tive to a deformation of the band structure provided the gap is 
open. On the other hand, it changes its sign as the Dirac mass 
changes its sign. Such a quantity is a topological charge. 
Hence an insulator phase is indexed by a set of four Chern 
numbers C^. A topological phase transition occurs when the 
sign of the Dirac mass A'] changes. 

It is instr uctive to make a reinterpretation of the Chern 
numbei^^®^. When the Hamiltonian is given in terms of 


the 2 X 2 Hamiltonian as in (2.9), or = r ' d with 


dx = rfhvYikx, dy = hvYfky, dz = A^, the Chern number 
is equivalent to the Pontryagin number. 




f dd A 


■d. 


(3.6) 


The Pontryagin number is a topological number which counts 
what times the vector d wraps a sphere. We use the polar 
coordinate of the d vector, d^Fidy = dz = 

cr{k), and we obtain 

^ y* d?k SijdicrdjO = - ^ J da, (3.7) 

which agrees with ( |3.5| ). The pseudospin texture forms a 
meron structure in the momentum space. A meron is a topo¬ 
logical structure which has a half integer Pontryagin number 
as shown in Fig|^ 



FIG. 6: Illustration of a meron structure in momentum space. A 
meron with the Pontryagin number (a) 1/2, which core spin is point¬ 
ing up direction, (b) — 1/2, which core spin is pointing down direc¬ 
tion. 


B. Classification of topological insulators 

We have defined four Chern numbers . Equivalently we 
may define the total Chern number C, the spin Chern num¬ 
ber t he val ley Chern numbe^^®^ and the spin-valley 

Chern numbeJ^^I^, 

C = (3.8) 

Cs = liC^ + Cf-C^-Cf), (3.9) 

= cf - cf + Cf - cf, (3.10) 

- cf - Cf + cf). (3.11) 

We make an important comment. The valley Chern number 
and the spin-valley Chern number are well defined only in the 
Dirac theory. Namely they are ill defined in the tight-binding 
model. Hence, we may call C and Cs the genuine Chern num¬ 
bers. 

Possible sets of genuine Chern numbers (C,Cs) are 
(0, 0), (2, 0), (0,1), (1, |) up to the sign ±. They are the 
trivial, quantum anomalous Hall (QAH), quantum spin Hall 
(QSH), spin-polarized quantum anomalous Hall (SQAH) in¬ 
sulators, respectively. Note that there are two-types of triv¬ 
ial band insu lators, which are quantum valley Hall (QVH) 
insulatoi^^^I^, and quantum spin- vallp^ Hall (QSVH) insula¬ 
tor with antiferromagnetic (AF) ordeW 

We comment on the relation between the Z 2 index and the 
spin Chern number. The spin Chern number Cg is identical to 
the Z 2 index by modulo 2 when there exists the time-reversal 
symmetr}!^. The spin Chern number is well defined even 
when there is no time-reversal symmetry, while the Z 2 index 
is well defined even when 5;^ is not a good quantum number. 


IV. TOPOLOGICAL EDGE 
A. Bulk-edge correspondence 

The most convenient way to determine if the system is topo¬ 
logical or trivial is to employ the bulk-edge correspondence. 
When there are two topological distinct phases, a topological 
phase transition must occur between them. It is generally ac¬ 
cepted that the band gap must close at the topological phase 
transition point since the topological number cannot change 
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cf 

cf 


cf 

c 

2Cs 

Cv 

2C,v 

QAH 

1/2 

1/2 

1/2 

1/2 

2 

0 

0 

0 

SQAH 

1/2 

1/2 

1/2 

-1/2 

1 

1 

1 

-1 

SQAH 

1/2 

1/2 

-1/2 

1/2 

1 

-1 

1 

1 

QVH 

1/2 

1/2 

-1/2 

-1/2 

0 

0 

2 

0 

SQAH 

1/2 

-1/2 

1/2 

1/2 

1 

1 

-1 

1 

QSH 

1/2 

-1/2 

1/2 

-1/2 

0 

2 

0 

0 

QSVH 

1/2 

-1/2 

-1/2 

1/2 

0 

0 

0 

2 

SQAH 

1/2 

-1/2 

-1/2 

-1/2 

-1 

1 

1 

1 

SQAH 

-1/2 

1/2 

1/2 

1/2 

1 

-1 

-1 

-1 

QSVH 

-1/2 

1/2 

1/2 

-1/2 

0 

0 

0 

-2 

QSH 

-1/2 

1/2 

-1/2 

1/2 

0 

-2 

0 

0 

SQAH 

-1/2 

1/2 

-1/2 

-1/2 

-1 

-1 

1 

-1 

QVH 

-1/2 

-1/2 

1/2 

1/2 

0 

0 

-2 

0 

SQAH 

-1/2 

-1/2 

1/2 

-1/2 

-1 

1 

-1 

-1 

SQAH 

-1/2 

-1/2 

-1/2 

1/2 

-1 

-1 

-1 

1 

QAH 

-1/2 

-1/2 

-1/2 

-1/2 

-2 

0 

0 

0 


TABLE II: Corresponding to the spin and valley degrees of freedom, 
there are 4 Chern numbers each of which takes zb|, Equiva¬ 
lently they are given by the Chern, spin Chern, valley Chern and 
spin-valley Chern numbers C, Cs, Cv and Csv They are indepen¬ 
dently controlled by the four parameters Aso, Av^, An and Xsx- 
Hence there are 16 states indexed by them. The genuin topological 
numbers are only C and Cg . 


its quantized value without gap closing. Note that the topolog¬ 
ical number is only defined in the gapped system and remains 
unchanged for any adiabatic process. 

To reveal the emergence of gapless modes at a phase tran¬ 
sition point, it is convenient to analyze the energy spectrum 
of a nanoribbon in a topological phase, because the boundary 
of the nanoribbon separates a topological state and the vac¬ 
uum whose topological numbers are zero. Indeed, we have 
pointed out the emergence of gapless edge modes in silicene: 
See Figj^bS). We may call it a topological edge when it sep¬ 
arates two topologically distinctive states. 


B. Helical edges and Chiral edges 

We analyze the edge modes of silicene in detail when sil¬ 
icene is in the QSH phase. The quantum numbers of the edge 
modes are not clear in Figj^bS) since they are degenerate. We 
may resolve the degeneracy by applying week electric field 
to a nanoribbon with zigzag edges as in Fig|7];a). It is an in¬ 
triguing feature of a zigzag nanoribbon that one of the edges 
is composed of A sites while the other of B sites. Hence 
the edge modes occurring in the A sites have higher energy 
than those in the B sites. This property enable us to iden¬ 
tify the edge modes occurring in the upper or lower edges. 
Second, the spin degeneracy is resolved due to the SO inter¬ 
action. We have shown up(down)-spin band in magenta(cyan) 
in Fig|^al). Third, the velocity of electrons within the edge 


(a1) 



FIG. 7: (al) Nanoribbon in QSH phase. There are helical edges. We 
have applied Ez to differentiate the upper (A sites) and lower (B 
sites) edges, which results in the splitting of the helical edges. (a2) 
Illustration of helical edge states. Up and down spins flow into the 
opposite directions, (bl) Nanoribbon in QAH phase. There are chiral 
edges. (b2) Illustration of chiral edge states. Up and down spins flow 
into the same directions. 


mode is identified by the slope of the band gap. In this way 
we know how the edge current fiows in Fig|7ja2). A promi¬ 
nent feature is the up and down spins flow into the opposite 
directions along each edge, implying that the edge current is a 
pure spin current. Such an edge is called a helical edge. 

We may also analyze a silicene nanoribbon in the QAH 
phase. The band structure is given by Fig|7Jbl), where we 
have applied week electric field Ez to resolve the degeneracy. 
The spin and the current direction are shown in both edges in 
Fig|7jb2). The edge current does not convey spins. Such an 
edge is called a chiral edge. 


C. Inner edges 

When a nanoribbon has only the valley Chern number and 
the spin-valley Chern number, no edge modes emerge because 
these numbers are not defined in the vacuum. Nevertheless 
these numbers are also topological numbers within the bulk. 

We may consider a junction separating two different topo¬ 
logical phases in a single honeycomb systenP^. We call such 
a junction an inner edge. In contrast we may call a real edge 
of a nanoribbon an outer edge. There is a crucial difference 
between gapless edge modes appearing along an inner edge 
and an outer edge. Any gapped state is indexed by a set of 
four topological numbers (C,C^^). Consequently, an 
inner edge state carries a gapless edge mode indexed by the 
difference (AC, AC^, AC^, ACg^) between the two adjacent 
gapped states. More precisely, we set AC = C^ — C^ and so 
on, when the topological insulator with (C^,Cj",C^,Cfi;,) is 
on the left-hand side of the one with (C^, C^, C^, Cf^). On 
the other hand, an outer edge state can carry a gapless edge 
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FIG. 8: (a) Silicene nanoribbon placed parallel to the x axis. When 
external field Ez is applied along the y axis, the region with Ez > 
Ecr becomes a QVH insulator. The inner edge between the QSH and 
QVH parts is helical. The outer edge below the QSH part is also 
helical. No gapless edge states appear along the outer edge above 
the QVH part, (b) The band structure of a silicene nanoribbon. Four 
gapless edge modes are found, which are assigned as in (a). 


mode only indexed by (C, C^) of the gapped state because the 
valley Chem numbers are ill defined in the vacuum. 

We illustrate a nanoribbon which contains the QSH and 
QVH phases in Fig[^ where there are two outer edges and one 
inner edge. As we have argued, the outer edge of the QSH part 
is helical, while the outer edge of the QVH part has no gapless 
edge modes. The nature of the inner edge is seen by analyzing 
the band structure in Figj^b). It contains only four nondegen¬ 
erate gapless states. On one hand, two solid lines correspond 
to the helical edge between the QSH part and the vacuum. 
On the other hand, two dotted lines correspond to the inner 
edge between the QSH and QVH parts, which is also helical. 
It should be noted that the inner edge modes cross the Fermi 
energy at the K and K' points. 



FIG. 9: (a) Y junctions of topological edge currents, which are made 
of the boundaries of three different topological insulators, (b) and (c) 
Topological circuits made of topological edges. We apply electric 
field to the center region of the QSH phase in (b). The region, to 
which electric field is applied, becomes a QVH phase. As a result, 
the connected edge states are separated as shown in (c). 


V. GAPLESS EDGE MODE IN DIRAC THEORY 

We co nstruct the Dirac theory of a gapless inner-edge 
stat^^l^ indexed by r] and s. It emerges along a curve where 
the Dirac mass vanishes, (x, y) = 0. 


A. Rectangular geometry 

We take the edge along the x axis. Zero modes emerge 
along the line determined by {y) = 0. We may set 
kx =constant due to the translational invariance along the x 
axis. The equation of motion Hk'iP = 0 reads 


D. Topological Kirchhoff law 


{-ihvYTydy + A^r^) '0 = 0. (5.1) 


We consider a configuration where three different topolog¬ 
ical insulators meet at one point: See Fig|^ In this configu¬ 
ration there are three edges forming a Y-junction. The condi¬ 
tion which edges can make a Y-j unction is the conservation of 
these topological numbers at the junction. This law is a remi¬ 
niscence of the Kirchhoff law, which dictates the conservation 
of currents at the junction of electronic circuits. We call it the 
topological Kirchhoff 

We present an interesting interpretation of the topological 
Kirchhoff law. We may regard each topological edge state as 
a world line of a particle carrying the four topological charges. 
The Y-j unction may be interpreted as a scattering process 
of these particles. In this scattering process, the topological 
charges conserve. We have shown that we can control the 
mass of Dirac cones with the spin and valley independently 
in silicene. Our findings will open a new way to topological 
spin-vaheytronics, where the spin and valley degrees of free¬ 
dom and the topology are fully manipulated. 


We set 0B = with ( = ±1, and seek the zero- 

energy solution. Here, 0^ is a two-component amplitude 
with up and down spins, 0^ = (0^,0^). Furthermore, 
setting -tpA {x,y) = (y), we obtain Hni^A {x,y) = 

Er](^'ipA y) together with Ey(^ = rjC^hvYkx, and 

{S,hv^dy - {y)) 4>A{y) = 0. (5.2) 


By solving this, the zero-energy solution is given by 


(a:, y) = Ce®*"* exp 




(5.3) 


and ipB {x^y) = (x^y), where C is the normalization 

constant. The sign of ( = ±1 is determined so as to make 
the wave function finite in the limit \y\ oo. This is a remi¬ 
niscence of the Jackiw-Rebbi mode^ presented for the chiral 
mode. The difference is the presence of the spin and valley 
indices in the wave function. 
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B. Circular geometry 


We consider a cylindrical symmetric domain with the radius 
ro at the origin of the xy plane. A phase transition occurs at 
r = ro, where /SJl{r) = 0. The equation of motion Hx'ip = 0 
reads 


+ (^idr = 

(^idr + V’s,^ “ = 

We can solve this for zero-energy states as 

-A f A^{r)dr' 

tlVF Jo 


C 




0 , 

0. (5.4) 


(5.5) 


and ipB (^ 5 ^) = K'^a (^ 5 ^), where C is the normalization 
constant and ( = ±1. The sign of ( is determined so as to 
make the wave function finite in the limit r ^ oc. 


C. Interface induced by electric field 

We apply an inhomogeneous electric field, 

A^(r) = Aso-Av^(r) (5.6) 

with Ay (r) = iEz (r) > Aso for r < vq. The region r < vq 
is a trivial insulator, while the region r > ro is a quantum spin- 
Hall insulator. The zero-energy helical edge current fiows 
along the circle r = vq. 


VI. TOPOLOGICAL QUANTUM FIELD-EFFECT 
TRANSISTOR 

We next calculate the cond uctance of a nanoribbon by using 
the Landauer formalisnl^^^. The conductance is quantized in 
silicene nanoribbons. Indeed, one channel has a quantized 
conductance e‘^/h. Accordingly, the conductance is obtained 
by counting the number of bands. We show the conductance 
in Fig[T^ When electric field is not applied, there are heli¬ 
cal edge states, which contribute to the conductance 2e^/h 
since up and down spin channels contribute to the conduc¬ 
tance. When the electric field exceeds the critical value 
Ecr, the edge states disappear since the nanoribbon becomes 
a trivial insulator, which results in zero conductance. This 
means the system acts as a transistor where the "on" state can 
be switched to the "off" state by applying electric field. This 
transistor is "quantum" since the conductance is quantized, 
which is highly contrasted with the ordinal transistor, where 
the conductance is not quantized. Furthermore the conduc¬ 
tance is topologically protected because the zero-energy edge 
state is topologically protected. Namely the conductance is 
robust against impurities due to its topological stability. Con¬ 
sequently we may call it a field-effect topological quantum 
transistoi^^. This is the most energy-saving device since it 
utilizes the minimum conductance. 


Band Structure 
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Conductance 
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o 8 8 88 
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FIG. 10: (Color online) Band structure, DOS and conductance of 
zigzag silicene nanoribbons for (a) the QSH insulator phase, (b) the 
metalic phase at the phase transition point, and (c) the trivial insulator 
phase. These phases are obtained by applying electric field Ez. The 
phase transition occurs at Ez = Ecr. The number of bands is 2W -h 2 
in the nanoribbon with width W. Here, the width is taken to be VF = 
31, and only a part of bands are shown. The band gap is degenerate 
(nongenerate) with respect to the up (red) and down (blue) spins at 
Ez = 0 {Ez > 0). Van Hove singularities emerge in the DOS at 
the points where the band dispersion is flat. The site-resolved DOS 
of the up-spin state at the outmost A and B sites of a nanoribbon 
are shown by red curves in the insets. There are flnite DOS for the 
zero-energy edge states in the QSH insulator. The conductance is 
quantized by unit of e^//i. 


We have calculated the density of states (DOS) p{E) and 
the conductance cr{E) of a nanoribbon as functions of the 
Fermi energy E, which is controlled by doping. We give the 
results at electric field Ez = 0, E^v and 2Ecv in Fig[^ A 
van Hove singularity occurs in the DOS at the point where the 
band dispersion is fiat. As E increases beyond the point, the 
Fermi level crosses a new band. A new channel opens and 
contributes to the conductance hy e^jh for each spin and val¬ 
ley. It is clearly observed that the edge channel connects the 
tips of the Dirac cones with the same spin at the K and K' 
points. 

We have also plotted the site-resolved DOS pi{E) of the 
up-spin states at the outmost A and B sites of a nanoribbon 
by red curves in the insets [Fig[T^. They represent degenerate 
zero-energy states dA Ez =0. As we have explained in Fig|7] 
the energy of the A and B sites become different for Ez ^ 0. 
It results in the downward (upward) shift of Pa{b){E) along 
the edge as Ez increases. They are separated completely, and 
zero-energy states disappear for Ez > F^cr- 
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VII. QUANTUM ANOMALOUS HALL EFFECTS 

Silicene has an additional interaction term, which we have 
so far neglected. It is the Rashba SO interaction given by 


(a) Band structure 


QAH 


(b) Berry curvature 

l^LLJuL 


(c) Spin 




E 




(cr X dij ^ Cjjj 

(7.1) -05 

Av / y -" ' - ■ N. 

0 

V J a 13 

-1 0 


■0.5 



in the tight-binding model. Here, Ar represents the Rashba SO 
coupling strength associated with next-nearest neighbor hop¬ 
ping, where /i^ = ±1 for the A (B) site, and dij = dij / \dij \ 
with the vector dij connecting two sites i and j in the same 
sublatticd^. The Dirac theory of the Rashba term is given by 

rjTzaXR {kyaj, - k^cjy ). (7.2) 

The Rashba term vanishes at the K and K' points. Hence 
the Rashba interaction is negligible as far as the low-energy 
physics near the Dirac points is concerned. 

There exists an exceptional case in which we cannot neglect 
the Rashba interaction. We have previously studied the inter¬ 
action term ( [2.13| ) affecting the Dirac mass. There is another 
type of interactions, which do not contribute to the Dirac mass 
but shift the bands into the opposite directions between up and 
down spins. The most important one is given by the exchange 
interaction. 


FIG. 11: (a) The energy spectrum of a QAH state. Gray curves are 
for a nanoribbon. Heavy green curves represent the energy spec¬ 
trum of the bulk, calculated independently. A gap opens in the bulk 
spectrum, where almost flat gapless modes appear at the edges of 
a nanoribbon. A red (blue) arrow indicates the spin direction away 
from the Fermi level, (b) (a) Berry curvature, (b) spin, and (c) band 
structure of a QAH state calculated based the Dirac theory. Gray 
curves represent the energy spectrum of a nanoribbon. Spins ro¬ 
tate by the Rashba interaction near the Fermi level, generating a 
skyrmion spin texture in the momentum space. It generates a non¬ 
trivial Berry curvature along the anticrossing circle. The integration 
of the Berry curvature gives the Chem number (7 = 2, since there 
are two skyrmions each of which yields (7 = 1. 


where k^, kp and kj run through ko = iuj, kx and ky, while 
G {uj, k) is the Matsubara Green’s function. 


M'^scl^Cis or Maz (7.3) 

is 


G{k) = 


1 

iuj-H {k) ’ 


( 8 . 2 ) 


in the tight-binding model or in the Dirac theory. We show the 
band structure of a nanoribbon with the exchange interaction 
in Fig[TTJa). We see that the two Dirac cones with the op¬ 
posite spins approach as |M| increases and touch each other 
at |M| = Aso- Then, we expect naively the level crossing to 
occur for |M| > Aso, that is, the two Dirac cones to penetrate 
into each other. Actually, the level crossing turns into the level 
anticrossing due to the spin mixing caused by the Rashba in¬ 
teraction Ar. The new phase is the QAH insulatoil^ with the 
Chern number 2. We show the Berry curvature in Fig[TT]^b), 
where it takes a large value where the bands almost touch the 
Fermi energy away from the K and K' points. The spin direc¬ 
tion is inverted and the spin forms a skyrmion texture, which 
results in the Chern number 1. The total Chern number is 2 
since there are two skyrmions at the K and K' points. 


VIII. SYMMETRY PROTECTED TOPOLOGICAL 
CHARGE 


with iuj referring to the Matsubara frequency (cc: real). This 
formula has a merit that it can be used even in the presence of 
interactionJ^^M^. 

We propose to define the following quantitieJ^, 



xTT[xiGdk^G-^Gdk,G-^Gdk^G-^i (8.3) 


for i = s,v,sv together with Xs = Xv = Vz, Xsv = 

cFzTjz. Provided the spin cfz and the pseudospin r]z are good 
quantum numbers, namely, [77, dz] = [77, r]z] = 0, we are 
able to prove that these are identical to ( |3.9| ), (|O0ll and 


by following the methocPl. We may use the formula (8.3 


even if symmetries are broken. They are no longer quantized 
when the symmetries are broken. Accordingly, we call them 
symmetry protected topological charges. 

In the presence of the Rashba interaction, the spin Chern 
number is explicitly calculated as 


The Chem number C is always quantized. However, the 
other ones Ci are not so when the relevant symmetry is bro¬ 
ken. Indeed, the spin Sz symmetry is broken by the Rashba 
interaction. We here analyze this problem. 

There exists an alternative expression for the Chern number 
in terms of the Green functiorPl. It reads 



xTr[(79/c^(7 ^Gdk^G ^Gdk^G ^], (8.1) 


[1 + {aXR/hvY^)‘^] ’ 

where Cs is the spin Chem number without the Rashba inter¬ 
action. This yields (3s = l — 5.9x10“^, where we have used 
= 5.5 X lO^m/s, a = 3.86A and Ar 2 = O.TmeV as sam¬ 
ple parameters of silicene. Surely Cg is not quantized, but the 
deviation is negligibly small. Furthermore, we can check that 
the edge mode remains practically gapless. The spin mixing 
can be neglected in practical purposes. 
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IX. TOPOLOGICAL SUPERCONDUCTOR 


Topological superconductors are superconductors which 
have nontrivial topological numberJ^. The bulk band spec¬ 
trum has a full gap due to the superconducting gap, while the 
edge states appear at the sample of edges as in the case of 
topological insulators. An intrig uing f eature is that these edge 
states support Majorana fermiong23ll3. 


A. Majorana fermion 

Majorana fermion is one of the hottest topics in the con¬ 
densed matter physics. Majorana fermion operator 7 satisfies 
the anti-commutation relation { 7 , 7 ^} = 1 and the Majorana 
condition 7 = 7 t, indicating that a Majorana fermion is a 
particle which is its own anti-particle. 

In the BCS theory, superconductor is expressed by the 
Bogoliubov-de Gennes (BdG) Hamiltonian, 

^ ^ E ( 4 . C-fe) ^fBdG (k) j , (9.1) 

based on the Nambu representation which is a 

combination of the particle and hole operators. The BdG 
Hamiltonian has a particle-hole symmetry, 

= —i^BdG, (9.2) 

with respect to the particle-hole operator E satisfying = 1 . 
The states with the energy E and —E are to be identified in 
the Nambu representation. The creation and annihilation op¬ 
erators of these two states satisfy 7 ^ = 7 -^;. Consequently, 
we obtain 7 q = 70 for a single zero-energy state (E = 0) if 
there is such one. As a result, a single zero-energy state with 
the particle-hole symmetry is a Majorana fermion. 


B. Topological superconductor in honeycomb system 

A topological superconductor is obtained from a QAH in¬ 
sulator by attaching an s-wave superconductor to iP^. In¬ 
deed, Cooper pairs are formed between up and down spins 
at the same site of the honeycomb system. The resultant BCS 
Hamiltonian reads H^cs = Hk + Hk' + i^sc with 

Hsc = E [^scc 7 ik) cfy {-k) + Asccf^^ {k) cf/ {-k)] 

r=A,B 

+ h.c. (9.3) 


this Hamiltonian. The BCS Hamiltonian is rewritten into the 
BdG Hamiltonian, 


^BdG = 


Hk ik) \ 

Hi i-k) ) ’ 


(9.4) 


by introducing the Nambu representation for the basis vector, 
i.e., «' = • 

Diagonalizing the BdG Hamiltonian, we obtain the energy 
spectrum. It consists of eight levels with the eigenvalues 


E^iik) = ±\l{nv,k)^ 



(9.5) 


with 

E’o ({Xso — olXv)^ + ^Ic) + P + Q(Asx) , 

(9.6) 

where a and P takes ±1. The gap closes (£’q= 0) at 

(Ah + Q^Asx)^ = (Aso — otXy)^ + (9.7) 

Though the original Hamiltonian is an 8 x 8 matrix, we may 
decompose it into 4 independent 2x2 Hamiltonians, 


77^’^ {k) 


PE^^^ hv^k- \ 


(9.8) 


corresponding to a, P=_^±l. This Hamiltonian reproduces 
the energy spectrum \9.5\ . We may interpret PE^'^ as the 
modified Dirac mass due to the BCS condensation. 

It is straightforward to calculate the Chem number of the 
superconducting honeycomb system H^^q. It is determined 
by the sign of the modified Dirac masJ^, 

C = \E sgn(/3K’^). (9.9) 

a,/3=±l 


The condition of the emergence of a topological supercon¬ 
ductivity is C 7 ^ 0. Note that it is zero when the time- 
reversal symmetry is present. In order to obtain a non-zero 
Chern number, Asx or Ah must be nonzero. It should be no¬ 
ticed that the spin Chern number is no longer defined due 
to the BCS condensation of the up-spin and down-spin elec¬ 
trons. The topological phase diagram is easily constructed in 
the (Aso, Ay, Asx, Ah, Asc) space. The phase boundaries are 
given by The Chern number is determined from ([9^). 


C. Majorana fermion in honeycomb system 


in the momentum representation, where Asc is the supercon¬ 
ducting gap. A finite gap present in a superconducting state 
allows us to evaluate the Chern number of the state to deter¬ 
mine whether it is topological. Alternatively we may examine 
the emergence of gapless edge modes by calculating the band 
structure of a nanoribbon with zigzag edge geometry based on 


We consider a disk region in a honeycomb sheet, as illus¬ 
trated in Fig[^ We may tune parameters Aso, Ay, Asx, Ah 
and Asc to become space-dependent so that the inner region 
has a different Chern number from the outer region. There 
appears gapless edge modes at the phase boundary, where 
’^(r) = 0. According to a general theorem, as we have 
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FIG. 12: Illustration of a Majorana zero-energy state between two 
topological phases. By applying electric field Ez locally, we may 
create a topological spot {C = 1) in a trivial superconductor {C = 
0). There appears a zero-energy Majorana state at the phase bound¬ 
ary. 

reviewed, a zero-energy state with the particle-hole symmetry 
is always a Majorana fermion. 


There are several ways to make ’^(rc) = 0, since there 
are four independent mass parameters Aso, Ay, Asx, Ah and 
one superconducting gap Asc- A simple way is to keep only 
one parameter and Asc- We consider the case where electric 
field is applied to a disk region of an antiferromagnetic topo¬ 
logical superconductor (Asx 7 ^ 0 ), as shown in Fig[T^ Very 
strong electric field can be applied experimentally by an STM 
probe. We assume electric field is strong enough so that 

\v{r) = £E,{r) = ±Aso + (9.10) 

The critical field is of the order of Ef = 0.1VA“^. The 
inner region of the circle have a nontrivial Chern number C = 
1 and becomes a topological superconductor. On the other 
hands, the outer region of the disk have C = 0 and remains to 
be the trivial superconductor. As a result, there emerges one 
Majorana fermion localized at the boundary of the circle. Its 
wave function is given by 
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